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Abstract 

We consider tan /^-enhanced quantum effects in the minimal super symmetric stan- 
^ . dard model (MSSM) including those from the Higgs sector. To this end, we match 

the MSSM to an effective two-Higgs doublet model (2HDM), assuming that all 
SUSY particles are heavy, and calculate the coefficients of the operators that van- 
ish or are suppressed in the MSSM at tree-level. Our result clarifies the depen- 
dence of the large-tan (3 resummation on the renormalization convention for tan /3, 
and provides analytic expressions for the Yukawa and trilinear Higgs interactions. 
The numerical effect is analyzed by means of a parameter scan, and we find that 
the Higgs-sector effects, where present, are typically larger than those from the 
"wrong-Higgs" Yukawa couplings in the 2HDM. 



1 Introduction 



It is well-known that loop effects in the minimal supersymmetric extension of the stan- 
dard model (MSSM) can become large when the ratio tan/5 = v u /v d of the two Higgs 
vacuum expectation values is sizeable. An important example is the Higgs coupling to 
bottom quarks, Hbb [1,2]. In the limit of heavy superpartner particle masses, the relevant 
couplings are [3] 



This includes a coupling dyt to the "wrong" Higgs doublet, which is forbidden by super- 
symmetry at tree level. The quantity Syj, is an ordinary one-loop effect and not enhanced 
in any way. However, due to the large ratio of vacuum expectation values (vevs) the 
bottom quark mass 



receives large tan /^-enhanced corrections to its tree- level value m b = y b v^ [4] . So does 
the Hbb coupling, when yj, is eliminated in favour of m& as is usually done. In general, 
while the loop effect is not large by itself, it provides large corrections to couplings or 
observables, which are suppressed by the small vacuum expectation value v d or which 
vanish at tree level. This type of relative enhancement is clearly a one- loop effect, which 
does not repeat itself in higher loop orders. 

Here we investigate other effects of this type from the Higgs sector, where loop effects 
can lead to a mixing of the H u and H d fields through two- and four-point interactions, 
which indirectly modify the Yukawa couplings to the physical Higgs bosons. We also 
address the issue of renormalization conventions for tan j3 and their stability against ra- 
diative corrections, as investigated in a numerical approach in [5,6]. The indirect effects 
from the Higgs sector have not been considered explicitly in the literature, although our 
results are implicit in one- loop calculations of MSSM Higgs decay rates [1,2]. However, 
these calculations are usually available only in numerical form. Having simple analytic 
expressions at hand is often useful to acquire a better understanding of parameter de- 
pendences, especially for the MSSM with its large parameter space. Higgs-sector effects 
from two-point interactions have been considered recently in [7], while their consequences 
for BB mixing are discussed in [8]. 

The theoretical issues can be exposed most clearly in the decoupling limit, assuming 
that all superpartner particles are heavy, of order Msusy, while the standard model 
(SM) particles together with all five physical Higgs bosons are light, at most 0(M EW ), 
the electroweak scale. In practice, the decoupling limit may be a good approximation for 
MgusY as small as a few hundred GeV. The assumed hierarchy implies a certain amount 
of fine-tuning, since the soft SUSY breaking parameters and the \i parameter must be 
O(Msusy), while rh\ d = |/i| 2 + rri^ d are 0(M EW ) to make the Higgs bosons light. The 



y D H d eQD - 5y D H^eQD + h.c. 

-> -(VbVd + Sy b v u ) b R b L - y b H d b R b L - 5y b H°*b R b L + h.c. 
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Higgs-mixing parameter b must even be much smaller than M| w to achieve large tan /3. 
The first fine-tuning is the well-known "little hierarchy problem" , here built in by hand 
by the assumption M SU sy ^> Af EW , but the second is an extra fine-tuning of b specific 
to the large-tan (3 scenario, which has received less attention up to now. To simplify the 
discussion we consider only the third family of fermions and assume CP conservation 
in the MSSM. Including the CKM matrix and CP-violating MSSM phases would be 
straightforward. Under these assumptions we match the MSSM to an effective general 
two-Higgs doublet model (2HDM), and determine the physical Higgs fields, the Yukawa 
couplings, and trilinear Higgs couplings from the one-loop corrected effective 2HDM. 

2 Matching the MSSM to the 2HDM 

The Higgs and Yukawa terms of the MSSM Lagrangian are given by 
£ = (D,H d Y(D^H d ) + (D^H u y(D»H u ) - m 2 d H d H d - m 2 u HlH u 

- b (H d eH u - fljeflj) - | (H d H u )(W u H d ) - l -{g 2 2 + gf) (H d H d - W U H U ) 2 

- y t H u eQU + y b H d eQD + y T H d eLE + h.c. (3) 

The conventions are as follows: the antisymmetric 2x2 matrix e is defined with e±2 = — 1. 
The hypercharge U(l)y and SU(2) couplings are g± and g 2 , respectively. The quark 
fields are Q a = (t La ,b La ) T , U a = (t* R ) a , D a = (b* R ) a , the leptons L a = (u La ,r La ) T , 
E a = {T R ) a , where a = 1, 2 denotes a Weyl spinor index. For the Higgs doublets we also 
use the notation $i = H d , $ 2 = = eH*, such that $i and $2 both have hypercharge 
-1/2. Our MSSM conventions follow [9] except for y b = -[y^Rosick, y T = -[z/r]Ro S iek, 
b = — [rra^ 2 ] Rosick and A t = — [A t ] Rosick . In this way our definitions of the Yukawa couplings 
and A parameters are consistent with the SPA convention [10]. 



2.1 Parameters and renormalization of the MSSM 

The relevant parameters of the MSSM Higgs sector are g±, g 2 , m 2 d , b, y^t and the two 
vevs v U)d . It appears to be contrary to the spirit of interpreting the MSSM as the high- 
scale theory and the 2HDM as the low-energy effective theory to introduce the vevs as 
MSSM parameters, since the vevs should be found a posteriori by minimizing the 2HDM 
effective potential. However, since we wish to use standard MSSM renormalization 
conventions, and moreover introduce tan (3 as one of the MSSM parameters, we substitute 

in the Lagrangian (3). In this way, after further fixing the renormalization convention, 
the parameter tan/5 = v u /v d is well-defined, contrary to the general 2HDM, where tan/3 
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is unphysical, since it can assume any value by performing a rotation in the ($i,$2) 
fields. In the MSSM such rotations are excluded by holomorphy of the superpotential. 
An equivalent set of parameters consists of 

e,M w ,M z , M A , tan f3,t u>d ,m bjt , (5) 

related to the previous set by 

e _ 9i92 

9i + 92 



M w = ^=v, M z = -j= s Jg 2 1 +g 2 2 v, 

M\ = mt + m d , tan (3 = — , 

v d 

t d = m 2 d v d + - (gl + g\) v d (vj - v 2 .) - bv u , 

t u = rh 2 u v u + - (g\ + gl) v u (v 2 u - v 2 d ) - bv d , (6) 

with v = ^v 2 + v\. The quark masses are given by = yb.tVd.u at tree level. For 
vanishing tadpole couplings t U)d) Ma is the mass of the CP-odd neutral Higgs boson. 

The parameters so far are bare parameters. We introduce renormalized parameters 
by substituting p — > po — > p + Ap, where nowpo (p) is the bare (renormalized) parameter 
and Ap the corresponding counterterm. The renormalization conventions are specified 
by imposing conditions on the parameters of the second set: e, Myy, Mz, Ma are de- 
fined by the conventional on-shell renormalization conditions for the electric charge and 
the particle pole masses, the tadpole couplings t u<d by the condition that the tadpoles 
(one-point functions) vanish. The renormalization condition for tan/3 and the Yukawa 
couplings, and the definition of beyond tree level will be discussed later. The coun- 
terterms Ap for the parameters of the first set are then defined by requiring that (6) 
holds for the bare as well as the renormalized parameters. Specifically, for tan/5 this 
implies 

tan/3 + Atp = [tan/3 = j-^r- = — r — = 1 + . . . . (7) 

[Vd\o v d + Av d v d v d V v u v d J 

Employing tan/3 = v u /v d , the one-loop relation 

At, = tan/3 (^-^) (8) 
\ v u v d J 

follows. Similarly, one obtains a set of equations relating the shifts of all parameters of 
the second set to those of the first. 
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We do not need to discuss the renormalization of superpartner masses and the A 
and \x terms, since these occur only in loops. In a one-loop calculation, it is sufficient to 
know these parameters at tree-level. Other SM parameters like the strong coupling or 
CKM matrix play no role here. For the subsequent discussion of tan /5-enhanced terms 
field renormalization in the MSSM can also be ignored. Thus, we shall not introduce 
renormalization factors Z u ^ for the SU(2)-doublet fields. The divergent parts of the 
vev counterterms Av u ^ are identical to those in Z Utd , since the theory in the broken 
phase can be rendered finite by the same counterterms as the unbroken theory before 
the shift (4). Thus, when we ignore field renormalization, we should formally consider 
Av U) d and consequently Atp as finite quantities, which indeed they are in the large- 
tan j3 approximation, whose values are nevertheless fixed by the imposed renormalization 
conventions. 



2.2 The effective 2HDM Lagrangian 

We integrate out the heavy SUSY particles and expand the MSSM effective action in 
local operators. The result is a general 2HDM, where the Higgs- Yukawa terms read 

Ahdm = Z ab {D^ a )\D^ b ) -m^t$ 6 



y(^i) 2 
A 



i^ 2 ) 2 - A 3 ($I$i) ($ 2 $ 2 ) - A 4 ($I$ 2 ) ($5$!) 



(9) 



y ($I$ 2 ) 2 + (A 6 $I$i + A 7 $|$ 2 ) $I$ 2 + h.C. 

+Y b $ 1 eQb + Y t ^ 2 Qt + Y T ^ x eLr - 5y b $ 2 eQb + SyMQt - 5y T <$> 2 eLT + h.c. 

neglecting operators with dimension higher than four. The last line equals 

Y b H d eQb - Y t H u eQt + Y T H d eLr - 5y b H^'eQb + by t H c d eQt - 5y T H^eLr + h.c. 

in terms of the H u d fields, so the Yukawa couplings take their standard form. Since 
we assumed CP conservation in the MSSM, the matrices Z ab , m ah (a,b = 1,2) and 
the couplings Aj are real. Note that the kinetic terms of the Higgs fields are not yet 
canonical, because we match 1PI Green functions. The kinetic terms of the fermions 
can be disregarded, since their renormalization is an ordinary one-loop correction to an 
unsuppressed tree term. At tree-level 
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Figure 1: Diagrams entering the matching of m? 2 and Z\ 2 - 



In the following the one-loop corrections to the couplings that vanish or are small at 
tree-level are of particular interest: SZ X2 = Zi 2 , Sm 2 2 , 5A 5>6i7 , and Sy^t,r- The result of 
the calculation is summarized in the next section. 

2.3 One- loop results 

All matching calculations are performed in the SU(2)xU(l)y symmetric "phase" of the 
MSSM. The couplings zu and m\ 2 follow from an expansion of the off-diagonal H u Hd 
two-point function in its external momentum. The relevant one-loop diagrams with 
heavy superpartner particles in the loop are shown in Figure 1. For the kinetic-mixing 
term we find 

z 12 = /I* ^3y t A* H 2 (m iR ,m$) + 3y b A* b H 2 (m~ bR , m§) + y T A* T H 2 (m fR , m L ) 

+ 3g 2 2 M 2 H 2 (\fi\,M 2 ) + g\M x H 2 {\^M X ) j (11) 

For the off-diagonal mass term, we write m\ 2 = b + 5b + [A6]h ea vy and obtain for Sb the 
expression (11) with the replacement H 2 (mi,m2) — > —l2(Trii,m 2 ), where 

7 2 (mi, m 2 ) = ( - - In ^- + 12(2?) J , 

l07T z \e v A J 

i 2 {x =1 + 1 , 12 

1 — X 



H 2 {m},m 2 ) = ^h 2 (x), 

16tt 2 2m? V ; 



1 - x 2 + 2x\nx 
2(X) = [l^xf ' ( ' 

and x = m\/m\. The kinetic-mixing correction z 12 is finite, as it should be, while the 
correction to m\ 2 is ultraviolet divergent and requires regularization (here dimensional 
regularization, space-time dimension d = 4 — 2e and renormalization scale v). The 1/e 
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pole is cancelled by the counterterm contribution [Aft] heavy to m\ 2 . The counterterm 
contribution is the difference between the MSSM and the 2HDM counterterm, which 
survives the matching relation, and cancels the singularity in the heavy-particle loops. 
In the following we shall drop the label "heavy" in the counterterm contributions. We 
remark that we have kept complex conjugation on the /i and A parameters, where re- 
quired, if they were complex. However, by our assumption of CP-conservation, \i and 
At t b jT are real. Here and below the masses that enter the arguments of the loop functions 
are the gaugino and sfermion mass parameters in the soft supersymmetry-breaking part 
of the MSSM Lagrangian and not the physical mass parameters. The physical masses re- 
ceive additional contributions at tree level proportional to the vevs v Uj d from electroweak 
symmetry breaking, which are of higher order in the Mew/Msusy expansion. 

For the four-point Higgs couplings As^j, which vanish at tree level, we obtain the 
ultraviolet-finite one-loop expressions 



5X 5 = -3(^AtyytK 2 (mQ,m iR )-3(fj,*At)Y b K 2 (mQ, 



m 



b R ) 



(fi*A* T )Y T K 2 (m L ,m 



TR , 



+ (/i*) 2 M 2 g\ K 2 (\/J I \,M 1 ) + 3(//) 2 M 2 g\ K 2 (\fi\,M 2 ) 
+ 2{^) 2 M X M 2 g\g\ K 3 (\^\, M lf M 2 ) , 



(14) 



5\ 6 = -3fi*A* b yl{j 2 (mQ,m~ bR ) + J 2 (m^,mg)} 
- n*A* yl {J 2 (m L , m fR ) + J 2 (m fR , m L )} 



+ WvtAt 



+ 3fx*y b A* b 



-a 2 - 9 ! 
9 2 3 



J 2 {rriQ,mi R ) + 
J 2 (mQ,m bR ) + 



4^ 2 



J 2 (m iR ,mQ) 



Mm- bR ,m, 



+ fi*y T A* T (--) {[-gl + gl] J 2 (m L ,m? R ) + [-2g 2 ] J 2 (m ffl ,m z )} 



-3^* \fi\ A* t y t K 2 (rriQ,m iR ) - 3// \A b \ A* b y b K 2 (mQ,m, 



-fi*\Ar\ A* T y T K 2 (m L ,m 



- 3// M 2 g\ L 2 (H, M 2 ) - // M x g{ L 2 (\ f i\,M 1 ) 



// (M 2 + M 1 )gigiL 3 (\ f j l \,M 2 ,M: 



i) , 



(15) 



S\ 7 = - 3fi*y t A; y 2 { J 2 (mg, m- tR ) + J 2 {m~ tR , m^)} 
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3n*y t A* t 



- 3n*y b A* b 



9 2 3 



J 2 (mQ,m iR ) + 
J 2 (mQ,m- bR ) + 



4^ 2 



J 2 (m iR ,mQ) 

2gl~ 



J 2 (m~ bR ,mi 



2 , 2 

-9 2 + 91 



J 2 (m z ,m ffl )+ -2g\ 



J 2 (m ffl ,m z 



- 3^* | A t | A* t y t K 2 {mQ, m- tR ) - 3^* \y h fi*\ A* b y b K 2 (mQ 



fj,*\y T /j,*\ A*y T K 2 (mi,mf 



TR - 



- 3fi* M 2 g\ L 2 (\fi\,M 2 ) - // M l9 { L 2 (\fi\, Ml) 

- n* (M 2 + MO g\g\ L 3 (\fi\,M 2 , Mi) . 



(16) 



The diagrams relevant for the computation of 5Xi are shown in Figure 2. The loop 
functions read 
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Figure 2: Diagrams entering the matching of 8X7. 



with x = 77i|/m^, ?/ = jrij/mj. 

The loop- induced Yukawa couplings to the "wrong" Higgs field are given 

g 

Syb = - ^*VbM 3 gp 3 (M 3 ,mQ,m iR ) - fi*y b y t A* t J 3 (\fi\,mQ,m iR ) 

3 1 

+ -fj,*y b M 2 g%J 3 {\{i\, M 2 ,mQ) + -//y&Mi glJ 3 {M u rrig, m~ bR ) 

+ -n* yb Mi glMHMt, m~ bR ) + -fjr yb Mi glMHMt, m Q ) , 



Sift = ^*ytM 3 glJ 3 (M 3 ,mQ,m iR ) + fi*y b y t A* b J 3 (\fi\,mQ,m bR ) 

3 2 

- -fj l * yt M 2 giJ 3 {\fJ l \,M 2 ,mQ) + -fi*y t M 1 gfJ 3 (M u mQ,m iR ] 

2 
3' 



-//•ytAfj J 3 (|/x|, Mi, m t ~J + -/i^tJlij ^J 3 (|/i|, Mx, m^) 



<tyr = ^*y r M 2 giJ 3 (\(i\, M 2 ,m L ) - n*y T M x g\j 3 {M x ,m~ L ,m f 



TR > 



s 



+ fi*y T M 1 g 2 1 J 3 (\fi\,M 1 ,m fR )- -ii*y T M x g\ J 3 (\[i\, M 1} m L ) , (23) 

where 

J 3 (m 1 ,m 2 ,m 3 ) 
J3(x,y) 

The expression for 5yb agrees with the well-known result for the tan /3-enhanced cor- 
rection to the bottom quark mass [3,4], when the bino contribution and terms of order 
M EW /Mgug Y are neglected. Our full results for Sijb/yb and Sy t /y t agree with the expres- 
sions for the correspondent quantities e + ey^ 2 and — [e' + e' Y yl] given in [7], which were 
written in terms of the rescaled A-parameters, Ab t t — > yb,tAb t t- Our result for 6X5 agrees 
with [8], where the bino contributions have been neglected. 

It is no surprise that the one-loop effects of interest are all proportional to the fi- 
parameter, since the off-diagonal two-point terms and the couplings are protected 
by a U(1)pq symmetry that is broken in the MSSM only by the fi and b term [4]. Note 
that the quantum correction to mf 2 is of 0(M| USY ), while the b (or rather mf 2 ) parameter 
was assumed to be of order M| w / tan/3. Thus, the large-tan (3 scenario requires more 
severe fine-tuning (by a factor of tan j3) than the one to maintain the "little hierarchy" 
Mew ^ Mbusy- It therefore appears less natural in this respect than the generic MSSM. 

3 Higgs couplings 

We are now ready to derive the one-loop corrected Higgs couplings, which are tan/3- 
enhanced relative to their tree-level expressions. For any quantity Q we can write its 
NLO MSSM expression as 

y = ^trcc H~ ^hcavy + flight) 

(25) 

where Ghc&vy is the contribution from the scale Msusy, in practice from loops contain- 
ing SUSY particles including the corresponding counterterms, and (/light the remaining 
contribution from the electroweak scale (standard model particles and Higgses). 

We will be interested only in effective interactions that vanish or are suppressed 
at tree-level by a factor of tan/3 relative to their natural size, and we shall compute 
them at one loop only at leading order in l/tan/3. In this case Gn g ht is suppressed 
relative to Qhcavy and can be set to zero. This can be understood from the fact that 
these interactions are protected by the PQ-symmetry mentioned above. This symmetry 
is broken strongly by the \x parameter that enters the high-scale contributions {?hcavy, 
but only weakly in the low-energy contributions by the small b (or rather m\ 2 i but in 
one-loop diagrams there is no distinction) parameter, b ~ M\/ tan (3. The low-energy 
contributions (?ii g ht all involve an insertion of b, or an explicit coupling to the small vev 
Vd ~ Mew/ tan j3, and are therefore negligible. On the other hand, the contribution 
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<?heavy from the superpartner particle loops is already encoded in the effective couplings 
of the 2HDM (9). We therefore conclude that the leading one-loop contributions to the 
tree-suppressed Higgs couplings in the MSSM are simply obtained by computing the 
physical Higgs fields with the general tree- level 2HDM Lagrangian (9). Subsequently, 
the specific expressions for the couplings z\2, m\ 2 , 8X5^7, and 8yb,t,r given in the previous 
section are inserted. 



3.1 Vacuum expectation value 

We first derive a relation between the shift (counterterm) of the vev v d and tan j3 and the 
MSSM parameter b, which feeds into all other Higgs couplings. The tadpole countert- 
erm is determined by requiring that the one-loop H d tadpole (1PI one-point function) 
vanishes: 

^H d = t<i + At d + [r^ d ] 0jhcavy + [rVj 0j i ight = 0. (26) 

The tadpole coupling t d is given by (6) and set to zero by our renormalization convention. 
The tadpole counterterm in terms of the shifts of the vevs is also obtained from (6) and 
given by 

At d = Am 2 d v d + m 2 Av d + ~{Ag 2 + Ag 2 )v d {v 2 d - v 2 u ) - Abv u - bAv u 

+ 4O? + gl> 2 d Av d - -(g 2 + g 2 2 )v 2 u Av d - -(g{ + g 2 2 )v d v u Av u . (27) 

The (unrenormalized) SUSY loop contribution to follows from the linear term in the 
2HDM Lagrangian (9) after shifting the neutral Higgs fields by the vevs, which leads to 

[r//Jo,hcavy = 8m 2 n v d - 8bv u + 8\xvl - Z8\ & v 2 d v u + (8\ 3 + <5A 4 + 8\ b )v d v 2 u - 8\ 7 v z u , (28) 

where 5mf 1: <5Ai_4 denote the SUSY loop contributions to the corresponding 2HDM 
parameters, which, however, will not be needed later on. The low-energy contributions 
can be neglected, as explained above, so the desired relation for Av d follows from At d + 

[r^Jcheavy = 0. 

This condition simplifies in the large-tan (3 limit. Applying v d <C v u , b <^ M| w , (27) 
reads 

At d w m 2 d Av d - Abv u - -{g 2 + g 2 2 )v 2 u Av d w M 2 A Av d - Abv (29) 
The second approximation uses 

2 _ 2(Mf - M 2 W ) 2 _ 2M 2 W 

6 = ^Misin(2/3), 

1 1 M 2 

m 2 = -M 2 A sm(2(3) tan/3 - -Af|cos(2/3) « M\ + 

10 



1 1 M 2 

m 2 u = -M\ sin(2/3) cot (3 + ~M\ cos(2/3) « — ^, 

^ = v cos /3, f u = v sin /3 ~ v, (30) 

which follow from inverting (6) for the renormalized parameters together with t u = = 
0. Similarly, 

[IVjo.heavy ~ Sbv ~ 5\ 7 V 3 . (31) 

Solving for Avj results in 



(5b + Ab)v + 5\ 7 v 3 
AVd = Ml • (32) 

We see that unless the counterterm Ab is fine-tuned the shift of the vev is of order 
^susy/^ew; which is very large relative to its assumed value of order M EW / tan/3. As 
discussed above, the shift of the vev is a finite quantity at leading order in large tan /3, 
since Ab cancels the divergences in 6b. Similar statements now apply to the shift of 
tan/3. Inserting (32) into (8), we find 



Av d l 2 n 5b + Ab + 5\ 7 v 2 

~Mj 



Atp ~ - tan 2 /3 — - = - tan 2 /3 — - 2 — + 0(tan (3). (33) 

"a 



The Higgs masses themselves receive no tan /3-enhanced one-loop corrections, so Ma 
corresponds to the physical mass of the CP-odd Higgs boson even at the one-loop level 
in our approximations. 



3.2 Renormalization of tan (3 

Contrary to the other parameters of the second set, there is no natural renormalization 
convention for the parameter tan f3. Gauge-independence and numerical stability of 
various schemes have been investigated in [5,6]. Quite generally, we expect problems 
with stability of observables against adding perturbative corrections, if the finite part of 
a counterterm is large, since this leads to a large shift between the tree-level and one-loop 
corrected values of parameters and large counterterm contributions to observables. 

From (33) we conclude that tan/3 receives large relative corrections, i.e. Atp/ tan/3 oc 
tan /3, unless 

5b + Ab + 5\ 7 v 2 = 0(1/ tan/3). (34) 

Some of the stability features at large tan /3 discussed in [5,6] can be understood on the 
basis of this equation. The DR scheme for tan (3 has Atp = by definition and fulfills 
(34) as does the DCPR definition [11,12], for which Atg/tan/3 oc 1. Indeed, these two 
schemes were found to exhibit good stability. On the other hand, the m3-scheme (DR 
for the b parameter), implies 



At ? _ + a [8b] aa + 5X 7 v 



2 



- — tan/3 —5 , (35) 

tan (3 Ml 
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Figure 3: Virtual Higgs boson exchange. 



which is very large, of order M| USY /M| W x tan (5. The Higgs mass and tadpole scheme [5] 
also violate (34), and are unstable. 

Good perturbative renormalization schemes should thus have 5b + Ab + 5\jv 2 = in 
the large-tan j3 limit. Since b ~ M| w / tan/? but Sb ~ M| USY , this implies fine-tuning of 
the counterterm for the b parameter as expected. 

3.3 Physical Higgs fields 

To determine the physical Higgs fields, we first bring the kinetic terms of the 2HDM (9) 
into canonical form by a field redefinition. This is not necessary to determine the Higgs 
masses, which receive no tan /3-enhanced corrections, since there are no tree-level sup- 
pressions, or observables that involve only highly virtual internal Higgs lines, as for 
instance in low-energy flavour physics. The latter statement can be illustrated with the 
help of Figure 3, representing part of a diagram with an internal Higgs propagator. The 
expression for this diagram using the Lagrangian (9) is i \{Zq 2 — m 2 ) -1 ]& a V^, where 
the propagator is a 2 x 2 matrix in Higgs doublet space. In low-energy scattering the 
momentum transfer satisfies q 2 <C m 2 , so the term involving the kinetic-mixing matrix 
Z a b in the 2HDM Lagrangian drops out at leading order in the low-energy expansion, and 
cannot affect physical observables. If one performs a field rotation to diagonalize Z a b, 
the rotation affects the vertices V a and mass matrix m 2 b in such a way, that once again 
the effect disappears at leading order in q 2 /m 2 . However, in this paper our aim is to 
calculate the couplings of on-shell Higgs bosons (q 2 ~ m 2 ) to quarks and the Higgs fields 
themselves, relevant to Higgs production and decay, so we need the fields corresponding 
to the physical states. 

The field redefinition that diagonalizes the kinetic terms (and includes the shift of 
the fields by the vevs (4)) is 



i 




No + i 



f)^ + I±^; 



u • 



0* 




1* 




1 — a 



-o* 

d i 



2 



Z12K 




(36) 
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where z a b is the one-loop correction to Z a b, and we used that the hermitian matrix Z a t 
is also symmetric due to CP-conservation. The diagonal coefficients zn, z 22 can be 
set to zero, since they are ordinary one-loop corrections to a non-vanishing tree term. 
The interesting terms are those that mix H d with the complex conjugate of H u . The 
arbitrary quantity a parameterizes a real field rotation in ($1, $2) space, which preserves 
the diagonal form of the kinetic term. We could set a = 0, but prefer to keep it to 
demonstrate explicitly the independence of physical quantities on a below. Note that we 
do not rotate the fields and then shift them by the vevs, since the vevs (and tan j3) have 
been defined as parameters of the MSSM Lagrangian before matching to the 2HDM. 

After substituting (36) into (9), we perform a unitary (in fact orthogonal, on account 
of CP-conservation) field rotation to diagonalize the Higgs mass matrix. The transfor- 
mation to the physical Higgs fields h°, H°, A , H ± , including the pseudo-Goldstone fields 
G°,G ± , is 

Im#°\ 1 / Sp + 5sp Cp + dcp 
Im H° d j y/2\- [cp + Sep] Sp + 5s p 

+ 




H+ \ I sp + 5sp Cp + Sep \ I G 
H?J ~ {-[cp + Scp] sp + 5sp \H+ 



ReH° u 
ReH° d 




(37) 



where 5sp, Sep, Ss a , 5c a parameterize the correction to the corresponding MSSM tree- 
level rotation, and we use the conventional notation = sin <fr, = cos <fi. We already 
incorporated here that the correction 5cp to the tree-level mixing matrix turns out to be 
the same for the CP-odd and the charged Higgs fields. The mixing angle a is given by 



tan 2a = tan 2(3. (38) 

Ml- Ml 17 



The correction terms 5sp, 5cp are of the size of an ordinary loop correction, and hence 
relevant only if the corresponding tree contribution is suppressed. This is the case for the 
off-diagonal elements, since cp oc 1/ tan/?. We therefore neglect the 5s p terms relative 
to sp ~ 1. For the off-diagonal correction we obtain 

1 + a 5b + Ab + 5X 7 v 2 
Sep = -— z 12 + ^ . (39) 

The second term vanishes in "good" renormalization schemes. 

In determining the correction to a, the cases Ma > Mz and Mz > Ma should be 
distinguished. In the following we discuss explicitly only the case Ma > Mz- The other 
case follows roughly (that is, up to some signs) from interchanging h° and H°. For large 
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tan/3 we have tan 2(3 ~ —2/ tan/3, so (38) implies that either c a or s a is small, unless 
Ma is very close to Mz to compensate the tan (3 suppression. For Ma > Mz the relevant 
case is s a small and negative and c a near 1. The correction to tan 2a is (valid for both 

cases M A > M z , and M A < M z ) 



+ mSi m 4 (4o) 

For M4 > Mz we have <5s Q = \5t 2a , 8c a « (relative to c Q ~ 1). The field redefinitions 
that lead to the mass eigenstates therefore read 

H ° = 72 ( Cah ° + [Sa + 5Sq] ^° + i[c/? + 5c/3] A ° + ' S/3G °) ' 

= (-[s Q + h° + c a #° + - i[ C/J + <5 C/3 ] G°) , 
#7 = spH" - [cp + Sep] G~, 

Ht = [cp + 5cp}H + + spG + . (41) 

We briefly comment on the case M A = Mz, since (40) appears to be singular for 
this parameter input. The masses used here are tree- level masses, since we keep loop 
corrections only when they are enhanced by tan/3 relative to the tree term. When 
Ma = Mz the CP-even Higgs bosons h°, H° are mass-degenerate and in fact physically 
indistinguishable. Eq. (40) is derived under the assumption that the difference between 
the diagonal elements of the neutral Higgs boson mass matrix is not small. Otherwise 
one must keep the one- loop correction to the diagonal elements in the calculation of St 2a , 
and the approximation discussed here cannot be used. Thus the resonant enhancement 
in (40) is a real effect as long as M\ — M§ is large compared to M| w x loop factor. 
However, if the widths of the Higgs bosons are so large that the h° and H° intermediate 
states are not separated by some measurement, nothing distinguishes the large tree- 
unsuppressed c a H° component from the smaller [s a + 5s a }h° component in (41) and 
there is no observable tan /3-enhancement in this measurement. In the following, we do 
not consider this degenerate case further. 



6t 



2a 



az 12 + 



Ml 



M\ + M 2 Z 
Ml-Ml 



Z12 



'Ml +Ml 



Ml 



Ml 



Sb + Ab + dXjv 1 



3.4 Yukawa couplings 

The Yukawa couplings and quark mass terms follow from inserting (41) into (9). The 
bottom mass term is of special interest. The relevant bilinear terms are 



C 



b,mass 



(Vb + 8y b + Ay b )(v d + Av d ) + 5y b (v u + Av u ] 



b R b L + h.c. 



(42) 
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The physical (pole) quark mass m& is obtained from this mass term plus low-energy self- 
energy contributions calculated in the 2HDM. As before we keep one-loop corrections 
only when they are enhanced relative to the tree expressions. This allows us to neglect 
the low-energy contributions (see [13,14] for the standard model electroweak loops) as 
well as Av u , and Sy b in (42). Ay b is the counterterm that cancels the divergences from 
superpartner loops in Syb- The finite part of Ayb defines the relation between the physical 
bottom quark mass and the renormalized Yukawa coupling yb- There is no reason to 
artificially define this finite part to be parametrically larger than the one of Syb, so we 
may drop Ayb as well. Introducing the definitions 

Sfjt $Vb 5y T 
e t = — , e b = — , e T = (43) 

Vt Vb Vt 

we obtain the mass terms 

£mass = -m t t R t L - m b b R b L - m T f R r L + h.c. (44) 

with 

m t = y t v u = ytvsp, 

m b = y b (v d + Av d ) + Sy b v u = y b vsp [- 1 + e 6 J . (45) 

ytanp + At^ J 

Note that we do not expand in Atp, since according to (33) it counts as a tan /3-enhanced 
loop correction relative to tan j3, just as e& tan (3 does relative to 1. No such enhancements 
are present for the top quark mass, which retains its tree expression. The expression for 
m T is obtained from the one for mb by replacing (yb, e&) — > (y T , e T ). 

Solving (45) for yb and eliminating y^t by mb,t, as it is conventionally done, we obtain 
the Higgs- Yukawa interactions 

£ Yuk = 1 [7-^-3 + e * + 5c P " ] ~^r Lz A {iirtLA ~ V2i R b L H^ 




+ - t R t L h° + " a ' ^ " ' 2 t R t L H c 

+ ii R t L G° -V2t R b L G + ^ 

^[tan/3 + At,] r _ v - 

V2v(l + e b [tanP + Atp])\ 



S a + Ss a — €b — ^-7rZ\2 t 



~ ~~ b R bi,hP + — b R bi,H 
sp sp 

+ e b + Sep + -^r^z 12 ) (ib R b L G° - V2b R t L G' 



Vtan/5 2 
+ h.c, (46) 
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where as before we neglect ordinary one-loop corrections relative to unsuppressed tree 
terms, and the r Yukawa couplings are given by the bottom ones after replacing all bot- 
tom parameters and fields by the corresponding expressions for tau. Several comments 
can be made. First, the explicit dependence on the arbitrary parameter a cancels with 
the implicit dependence in (39), (40). Second, our result is independent on the renor- 
malization convention for tan (3. A change of renormalization prescription is effected by 
a change in the finite part of Atp, or (33) Avj, or (32) Ab. From (39), (40), we therefore 
find that 

Sep 



[tan (3 + At 



^— -roup ,[s a + 5s a ] ,sp,c a (47) 

are invariant under a change of renormalization prescription that contains only tan/3- 
enhanced one-loop terms. 

Eq. (46) shows that a resummation of tan /5-enhanced terms is necessary only in the 
overall factor multiplying the bottom (and tau) Yukawa couplings [3], where the need 
for resummation arises from using the quark (lepton) mass as a parameter rather than 
the MSSM Yukawa coupling. All of the Yukawa couplings which are tan /3-suppressed at 
tree-level receive unsuppressed loop corrections, which are incorporated in (46) through 
the "wrong-Higgs" Yukawa couplings e t>biT , and the quantities 5cp, 5s a , z 12 from the one- 
loop effects in the Higgs sector. The loop-induced couplings may well exceed the tree 
terms, if tan j3 is large. 

To compare our result with previous results in the literature, see Eq. (21) in [3], we 
consider the bottom quark coupling to the lightest Higgs boson h°, which is the product 
of two factors 

m b [tan f3 + Atp] ^ s a + 5s a - e b - ^z 12 



\f2v (1 + e b [tan/3 + Atp]) sp 

With respect to the first factor, note that we do not discuss logarithmic one-loop effects, 
and hence do not distinguish the physical quark mass m b from the running MS quark 
mass fn b (Q) used in [3]. In reality, renormalization group evolution from m b to the 
electroweak scale is an important effect that should be taken into account as done in [3]. 
The main difference in the first factor constitutes the presence of the (finite) counterterm 
Atp, which must not be neglected in general, and which renders the resummed result 
manifestly scheme-independent. This clarifies that the previously known result where 
Atp = 0, is only valid in "good" renormalization schemes (34) where the one-loop shift 
of tan j3 is not itself tan /3-enhanced, as also observed recently in the related context 
of the muon Yukawa coupling [15]. The resummation formula from [3] without the 
counterterm contribution should not be used in other schemes, for which Atp is an 
important correction in the first factor. 

Using sp ~ 1 in the one-loop corrections, and (33), (40), the second factor can be 
rewritten as 

Tp~ eb+ Ml -Ml 212 ~ M\- M 2 Z ~ Ml -Ml Ml ' ^ ' 
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The first two terms reproduce the corresponding result in [3], but the Higgs-sector effects 
embodied in the last three terms have not been included there. Higgs-sector effects 
related to kinetic mixing have recently been computed in [7] . The DR scheme for tan j3 
is assumed there, hence the last term involving 5b + Ab + SX^v 2 oc Aig = vanishes. We 
find several discrepancies in the comparison with [7] . The third term proportional to z±2 
would be consistent with the result of [7] if z i2 = e GP , but our result for z 12 in (11) gives 
Z12 = — cgp- We also find that the corrections to the iitRH coupling are not included 
in [7]. The most important difference is that the fourth term involving 6X7 is missed, 
since it is assumed that Higgs-sector effects arise only from kinetic mixing, and not (also) 
from the Higgs self-couplings. The large quantum correction to the tree-level suppressed 
charged Higgs tRbiH + coupling is present in the discussion of the b — > 57 transition in 
[16,17]. These papers account for the gluino contribution to e t . The Higgs-sector terms 
involving 5c@ and Z12 are again not considered. 

More precisely, in the expression for the third term proportional to z\ 2 Ref. [7] has 
M? i0 /(M H0 — M^q), where we have M\j{M\ — M§). The two expressions are consistent 
with each other, since in the one-loop correction we may use the large-tan (3 limits of the 
tree-level mass relations, which give M h o = Mz, and M H o = Ma- In reality, the Higgs 
masses must receive large quantum corrections to satisfy experimental limits, which are 
not included in our equations, since they are not tan /3-enhanced. Since the apparent 
singularity for M4 = Mz arises from the degeneracy of the two CP-even Higgs bosons, as 
discussed above, it is plausible that Ml /(M H0 — M% ) provides a better approximation, 
even though formally the difference to M\j(M\ — M|) is a higher-order effect. 

In Section 4 we perform some representative numerical estimates of the new Higgs- 
sector corrections to the Yukawa couplings and compare them to the previously known 
term e&. 

3.5 Higgs self-couplings 

The Higgs self-couplings are obtained in the same straightforward fashion. Here we focus 
on the trilinear couplings, since the quartic self-interactions with coefficients induced by 
quantum effects will be too small to be measured in any foreseeable experiment. The 
detection of tan/?- or loop-suppressed trilinear couplings provides a challenge even to 
the ILC (see, e.g., [18]). 

The Lagrangian for the trilinear couplings 

£-HHH = £-HHH, large + &HHH, small 

(50) 

consists of a piece Chhh, large involving interactions that do not vanish at tree level when 
tan/5 — > 00, and another piece Chhh, small with interactions that do. The tree-level 
couplings in the MSSM are standard and will not be repeated here. We find a compact 
expression for the tan /3-enhanced one-loop corrections relative to the small tree-level 
interactions in Chhh, small, applicable to the case Ma > Mz: 

= -4 {(-l)(3Mi + 2Ml) C H h° 2 H° 
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+ (M 2 A C H - 5\ - 5X 7 ) (H° 2 + A° Z )H C 



+ 2 (Ml-2M^)C H -(S\ 6 + S\ 7 ) 



H°H + H~ 



+ 2 (Ml- M 2 Z ) C H h°A°G° - 25\ 5 H°A°G° - Ml C H H°(G° 2 + 2 G + G~) 
+ 2 (M 2 A - M 2 Z + M*y) C H h°(H~G + + , (51) 



where 



c 



1 M 



H 



M 2 A -M 2 Z \ v 2 



8b + Ab + 8X 7 v 2 z 12 



Ml 2 



+ SX 7 ). (52) 



4 Numerical estimates 



In this section we present numerical estimates of the Higgs-sector corrections by evalu- 
ating them for degenerate SUSY masses, and by performing a random parameter space 
scan. We are especially interested in comparing the new Higgs-sector effects to the known 
one-loop induced "wrong-Higgs" Yukawa couplings. 

4.1 Degenerate SUSY masses 

We begin the analysis of the size of the new terms in the Yukawa couplings relative to 
€b,t by assuming that all SUSY parameters are equal to a common mass scale Msusy 
and by allowing for arbitrary signs in the /i and A parameters, but not the gaugino 
masses. In the following we give analytic results for e b t jT and the quantities z 12 and 
<5A5_7, which determine the Higgs-sector corrections to the Yukawa couplings and the 
Higgs self-couplings, and numerical results to judge the importance of individual terms. 
To make contact with the literature, we introduce the rescaled A parameters A by 
Ai = UiAj,. In units of sgn(yU,)/(967r 2 ), we find 

/ 1 11 \ 

e b = 3y 2 sgn(A) + 32na s - 187ra cm — + 



/ t s 3.02 
9.65 + sgn(A f ) — — , 



( 1 5 

e t = - 3y 2 sgn(A b ) - 32rra s + 187ra em — 



si 27c 2 



2 14 

-9.77 - sgn(vV 



[50 c p ] 



2 



e r = -18vra em (\ - J-] = -1.80 (53) 
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for the "wrong-Higgs" Yukawa couplings, and 

Z12 = 3yt sgn(A t ) + 3y\ sgn(A b ) + y 2 T sgn(A r ) + 47ra em ( 4" + 4" 



3.02 ,~ , 2.14 ,~ , 0.26 

T 



1.45 + sgn (j 4«) -3- + sgB(A) ^ + sgB(^) ^ 



<5A 5 = sgn(/i) [ - 3y? - 3y\ -y* + 16vr 2 a^ m [ — + \ + 



2 



C W S W , 



sgn(/i) 0.712 



3.04 1.59 



•s 



1 



[50 cpY 



'/3 

5A 6 = 15 j/j sgn(A 6 ) - 3% 4 sgn(A t ) + 5y 4 sgn(A T ) 

37T« r 



'cm 



+ 327r 2 aL I — + — + 



3y t sgn(At) + 3y fe sgn(A 6 ) + y T sgn(A T 



3 1 2 



em I 4 1 4 1 2 2 

°iu Si) °tu L 'w. 



, 3.04 1.29 \ ,~ , / 0.91 7.62 



sgn(A T 



0.11 0.34 



[50 C/3 ] 2 [50 c Y 



5X 7 = 15y?sgn(A t ) - 3y^sgn(A b ) - y 4 sgn(A T ) 

+ 3 2 Qc 2 m (~%t 2 S g n (^) + %& sgn(A) + yl sgn(I 7 

9 o / 3 1 2 \ 
+ 32vr 2 a 2 m — + — + 



cm l 4 1 4 1 2 r 2 

01 ) '-'0/1 ^H) ^1! 



, 15.2 1.29 \ /7 , / 0.91 1.52 



,~ . / 0.11 0.07 \ 

+sgnW " p^J (54) 

for the Higgs-sector couplings. Our results for 5A 5 „ 7 agree with the special case given 
in [20], which assumed all squark masses equal to M SU sy and neglected the gaugino 
and slepton contributions. The numerical values in (54) have been evaluated with SM 
parameter input as follows: for the gauge boson masses we choose Mjy = 80.398 GeV 
and Mz = 91.188 GeV, the gauge couplings are fixed to a em = a^(Mz) = 1/127.9 and 
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Figure 4: Dependence of Higgs-sector couplings, normalized to e^, on tan/5. For the 
ratios Zvijt^ SXej/ei, the plots correspond to the sign choices {sgn(A t ), sgn(y4&)} = 
{+>+}> {->+}> {+>-}> {->-}• In tne P lots for <5-W e f> we adopt {sgn(/i), sgn(A)} = 
{+,+}, {-,+}, {+,-}, {-,-}• 



a s = a^ LS (Mz) = gf/(47r) = 0.118. The cosine of the weak mixing angle is given by 
c w = Mw/Mz- For the top Yukawa coupling we use y t = m t /(vsp) ~ 1.003/s^, where 
m t = 171.7 GeV and v = s w M w / y/2na em ~ 171.2 GeV. For the down-type and lepton 
Yukawa couplings we here use the tree-level relations = rrib/(vcp) ~ 0.0169/C/5 and 
y T = m T /{vcp) ~ 0.0102/C/3, where m b = 2.89 GeV and m T = 1746.24 MeV. All fermion 
masses correspond to the renormalized masses in the MS scheme at fi — Mz [19] rather 
than the pole masses. 

Since e& is nearly independent of tan/3 for tan/3 ^> 1, we show the tan/3 dependence 
of the ratios z^/ei, and 5A 5 „ 7 /e6 from (53), (54) in Figure 4 for different choices of the 
signs of the A-parameters. The slepton terms give a negligible contribution except for 
Z\2 at very large values of tan/3, so we always set the sign of A T to +1. For tan/3 
smaller than 20 - 30 the value of z^, 6X5^7 is determined by the A t terms and the 
gaugino contributions. For SX7, in particular, which feeds into the Yukawa couplings, 
the terms proportional to A t give by far the largest contributions, while the effect of 
At, becomes relevant for very large tan (3 values, reaching up to a 20% contribution. On 
the other hand we observe from the plots that the Aj, terms and a change in the sign 
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of Ab has a large effect on 5\q, which is driven by the sgn(Ab)/cp term in (54) with its 
sizeable coefficient. In general, under the assumption of a common mass scale for all 
SUSY parameters, we thus conclude that the quantities z u , 5\$ t Q and particularly 5Xj 
are comparable in size with the "wrong-Higgs" couplings e^t even though they do not 
receive gluino contributions. Note that \5\f\ is larger than je^l for all choices of signs of 
the A tl b- 



4.2 Estimates using random parameter-space sampling 

We now investigate in more detail the size of the Higgs one-loop corrections to the 
Yukawa interactions relative to e b) t for different SUSY parameter input. We consider 
good renormalization schemes for tan/3, i.e. schemes satisfying 5b + Ab + 5Xrv 2 = in 
the large-tan (3 limit, in which case the two relevant corrections from (46) are 

h ~z h ~ M| * 12 " 2vHXl (55) 

h=Z 12 , h 2 = M 2 a _ M 2 ■ ^ 

In terms of the quantities hi 2, the one-loop improved couplings of the Higgs fields to 
the bottom and top quark are rewritten as, omitting the global factors outside the curly 
brackets in Eq. (46), 

b R b L h° : — + 5s a - e b - ^—^z 12 = — - e b + h 2 , 

Sp 2 Sp 

■T , ^yO 1 1 + a 1 

%b R b L G : + e b + dcp H —z 12 = + e b , 

tan p 2 tan p 

At) 1 e 1 — a 1 , 

it R t L A : + e t + dcp ~ z i2 = 7 75 + e t - h i > 

tan p 2 tan p 

i R t L H° : ^ + Ss a - e t + l -^z l2 = fSL - £t + hl + h 2 . (56) 

The function Ch, which appears in the one- loop corrections to the small tree- level tri- 
linear Higgs couplings (51), also has a simple expression in terms of h 2 , given by 

We perform a random parameter sampling for the SUSY parameters in the following 
ranges: 500 GeV < rhi < 5 TeV, with being the squark, slepton or gaugino masses, 
and 500 GeV< \A t \, \Ab\, \A T \ < 5 TeV. The different sign assignments for the pa- 
rameters fi, A t and A b are explored separately, but we fix sgn(A T ) = +1. Since the 
Ma dependence of h 2 comes in only through the prefactor 1/(M\ — Mf), we shall fix 
it to the reference value = 200 GeV. Here contrary to the analytic expressions in 
the previous subsection, we correctly include the tan /3-enhanced loop correction to the 
down-type quark and lepton Yukawa couplings, i. e. yb, T = ^&,t/ ( vc p) x 1/(1 + e f>,r tan j3) 
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Figure 5: The corrections hi and h 2 versus e t and e& respectively, for = 200 GeV and 
positive fi. The plots in the first row have sgn(A t ) = +1 and those in the second row 
have sgn(A t ) = — 1 (sgn(Ab) = +1 in all plots). 



is used. We remark that not all points allowed in the parameter scan are physical, since 
vacuum stability requirements do not allow large A terms relative to the scale of the 
SUSY particle masses. 

Let us first study the case of positive /1, where e b is preferentially positive, so that 
Ub is reduced by the one-loop correction. In Figure 5 we show hi (left) and hi (right) 
versus e t and e&, respectively, for both signs of A t (positive in the upper plots, negative 
in the lower ones). For positive \i, the value of tan/3 does not significantly change the 
shape of the scatter plot, and hence has been fixed to tan (3 = 35. The typical size of the 
effective couplings is a few percent. Values of \h 2 \ range up to 0.1, which implies that 
the one-loop induced coupling may exceed the tan /3-suppressed tree coupling for values 
of tan/3 as small as 10. In comparison, hi, which originates only from kinetic mixing, 
tends to be significantly smaller. Thus, the Higgs sector correction h 2 competes in size 
with the "wrong-Higgs" Yukawa coupling e& even for large values of Ma, and is in fact 
the dominant one-loop correction to the bbhP Yukawa interaction. Some general features 
that can be observed in Figure 5, or extracted from the analytic expressions are: 

i/ In the assumed range of SUSY parameters, the sign of e& (et) is correlated (anti- 
correlated) with the sign of \i, since the gluino contribution (first term in (21)) 
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dominates and the function J3 is always negative. The term proportional to A t 
can overcome the gluino contribution and change the sign of e& only for negative A t 
in the parameter-space region where |v4j|//z 2 is much larger than 1/M 3 . Similarly, 
large negative Ab can make e t positive, when yl\Ab\/ [i 2 ^> 1/M3. 

Since cannot reach large negative values, there is no significant enhancement of 
the bottom Yukawa coupling from the large-tan (3 resummation. The numerical 
effect of the terms proportional to y 2 Ab is therefore small compared to those pro- 
portional to y 2 A t . We verified that the scatter plots for negative Ab are similar to 
the ones with positive Ab shown in Figure 5. 

Larger values of \e t \ (|e&|) tend to be correlated with larger values of \h\\ (|^2|)- 
However, there is no strict relation, since the dominant contribution to the "wrong- 
Higgs" Yukawa couplings is proportional to the gluino mass, while those to the 
Higgs-sector couplings involve the A terms. 

The largest contribution to h\ is given by the term proportional to A t in (11). 
However, the sign of hi is not determined uniquely by the sign of A t due to the 
gaugino term proportional to g 2 , which is always positive (the function H 2 is posi- 
tive) and can become comparable to the A t term in some regions of the parameter 
space. 

In most of the explored parameter space the sign of h 2 is opposite to the sign of 
At. Since the 6X7 term in the definition (55) of h 2 is the larger of the two, we 
have sgn(/i 2 ) = — sgn(5Ar). The relation sgn(5Ay) = sgn(A t ) can be understood as 
follows: neglecting the smaller terms involving gauge couplings, Eq. (16) for 6X7 
contains the terms 

-3/iyfAt {J 2 (mg,m tfl ) + J 2 (m iR ,m§)} - 3fiAfyfK 2 {mQ,m iR ) (58) 

proportional to powers of A t . The first term has the same sign as A t , since J 2 
is always negative, while the second term has opposite sign, since K 2 is positive. 
Defining the ratio 



r(x) = _U^m R ) + Um R ,m QJ 



>-R ' 



with x = m 2 R /rriQ, we find that r(x) is monotonically increasing, satisfying 1.6 < 
r(x) < 169 when 1/100 < x < 100 as allowed by our parameter-space sampling. 
The third term in (58) can compete with the first only if (At/run) 2 > r(x). For 
x = 1 this relation requires large A terms, (A t /rriQ) 2 > 6, which are disfavoured 
by vacuum stability arguments. We also note that the largest gaugino term in 6X7 
(the one proportional to g 2 ) gives a positive contribution because the function L 2 
is negative. 
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Figure 6: The correction h 2 versus e b , for Ma = 200 GeV and negative The plots in 
the first row correspond to {sgn(y4 4 ), sgn(A b )} = {+1,+1} and { — 1,+1}, while those 
in the second row have {sgn(A t ), sgn(A b )} = {+1, —1}, { — 1, — 1}- Black, lightgray and 
darkgray points (blue, yellow and red in colored plots) correspond to tan/? = 10, 35 and 
60, respectively. 



Let us now turn to the case sgn(/i) = — 1. The most important difference is that the 
negative sign of the \x parameter flips the sign of e b to negative values (and the one of e t to 
positive values), which leads to a strong increase of the bottom and r Yukawa coupling, 
when e b cancels the 1/ tan/3 term in (45). The validity of a perturbative expansion is in 
doubt when the Yukawa couplings become too large. In the parameter space sampling 
for negative /i, we therefore only keep points that satisfy < y biT < 2. 

The dependence of the Higgs sector correction h 2 versus e b , for three representative 
values of tan/3 = 10,35,60, Ma = 200 GeV, and all other SUSY parameters scanned 
randomly in the above intervals is shown Figure 6. The perturbativity cut on y b has 
a strong effect, since it eliminates points with large negative values e b given a value of 
tan/3. The value of e b at which this happens can be estimated by e b tan(3 « —0.5, which 
gives — e b ~ 0.05,0.014,0.008 for tan/? = 10,35,60, respectively, in agreement with the 
figures. At these points we observe a rapid increase of 1^1 driven by the term 

-3/i 3 ytA b K 2 (mQ,m~ bR ) (60) 
in (16), which gives a contribution to h 2 ~ — SX-j with sign opposite to A b . Point 
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Figure 7: The correction hi versus e t , for Ma = 200 GeV and negative /i. Plots in the 
first row have (sgn(y4 t ), sgn(A b )} = {+1,+1} and { — 1,+1}, and {sgn(A t ), sgn(A b )} = 
{+1,-1}, { — 1,-1} in the second row. Black, lightgray and darkgray points (blue, 
yellow and red in colored plots) correspond to tan/5 = 10, 35 and 60, respectively 



i/ discussed above also applies to the case of negative /i taking into the account the 
reversed signs of e^. Finally, Figure 7 shows the the Higgs correction hi versus et for 
tan p = 10, 35, 60. Compared to the case of positive fi, the sign of A b has now a relevant 
effect, especially for larger values of tan/3, due to the contribution of the A b term in z 12 . 
However, the growth with y b is less pronounced than in case of h 2 , since hi does not 
contain terms proportional to y b . 

5 Summary 

This paper has been motivated by previous work [3] that systematically investigated 
the resummation of SUSY QCD large-tan P effects in the decoupling limit, where the 
standard model particles and Higgs scalars remain light, but which did not consider 
electroweak effects. We performed a complete one-loop matching of the MSSM to a two- 
Higgs doublet model for all couplings that are absent or suppressed in the MSSM with 
large tan P, keeping all contributions that are tan /3-enhanced relative to the suppressed 
tree terms. This includes complete expressions for the kinetic mixing and Higgs self 
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couplings 6X5-7, which have not been given previously. 

Our result confirms, as expected, that a resummation of large-tan (3 effects to ar- 
bitrary loop order is necessary only for the bottom and tau Yukawa couplings, if the 
bottom and tau mass are used as parameters of the MSSM. It clarifies a point left open 
in [3], namely how the renormalization of tan (3 affects this resummation. We find that 
the standard expression for the resummed Yukawa coupling is valid in renormalization 
schemes where the counterterms (shifts) of tan f3 and the vacuum expectation value Vd 
do not receive large-tan (3 contributions. In other schemes, a finite counterterm must be 
explicitly included in the resummation formula. 

Besides the known "wrong-Higgs" Yukawa couplings, all other Higgs-sector effects 
that feed into the Yukawa couplings by modifying the definition of the physical Higgs 
fields at one loop can be parameterized in terms of two couplings hi, h 2 dominated by 
kinetic mixing and the Higgs self-coupling 6X7, respectively. The same hi also enters the 
effective trilinear Higgs couplings. Our numerical study suggests that, where present, 
the Higgs-sector effect hi is more important than the "wrong-Higgs" Yukawa couplings, 
and we identified the dominant SUSY parameter dependences. Our result extends and 
corrects a previous result [7] that included the kinetic-mixing effect, but not the one 
from the one-loop induced Higgs couplings. 

Although obtained in the decoupling limit, the effective one- loop couplings derived in 
this paper should be useful to obtain simple analytic estimates of the leading quantum 
corrections to those Yukawa and Higgs trilinear interactions, that are suppressed at 
tree-level, such as the bbh° coupling. 
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